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ABSTRACT 

We show that compactly supported wavelets in L 2 (M) of scale N may be effectively parameterized with a finite set 
of spin vectors in C , and conversely that every set of spin vectors corresponds to a wavelet. The characterization is 
given in terms of irreducible representations of orthogonality relations defined from multiresolution wavelet filters. 
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1. INTRODUCTION 
r ~^\ Let L 2 (R) be the Hilbert space of all L 2 -functions on R. Let ip G L 2 (R), and set 

[JLr ip n ,k (x) ■= 2%^(2 n x- k) for x G R, and n,k G Z. (1) 
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We say that tjj is a wavelet (in the strict sense) if 

{ipn.k ;n,ke 1} (2) 

constitutes an orthonormal basis in L 2 (R); and we say that i/j is a wavelet in the frame sense (tight frame) if 

ll/ll^(M)= E K^l/>| 2 (3) 
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V.Q holds for all / G £ (R), where ( • | • ) is the usual L (R)-inner product, i.e. 
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{ipn.k I /) = / 1pn,k (x)f (x) dx. (4) 

Jm. 

It is knowrfl'i that a given wavelet ip in the sense of frames is a (strict) wavelet if and only if HV'IIl^r) = Either 
way, the numbers in (0) are the wavelet coefficients. 

We shall have occasion to consider scaling on R other than the dyadic one, say x h- > Nx where N G N, N > 2. 
Then the analogue of (Q) is 

ip ntk (x) :=N%ip(N n x-k), xeR,n,keZ. (5) 

However, in that case, it is generically not enough to consider only one ip in L 2 (R): If the wavelet is constructed 
from an 7V-subband wavelet filter as in Ref. |^, then we will be able to construct ip^ l \ip^ , . . . , i/'' Ar_1 - ) in L 2 (R) such 
that the functions in (|fj) have the basis property, either in the strict sense, or in the sense of frames. In that case, 
the system 

{ip® k ; 1 < i < N, «,fcez) (6) 

will constitute an orthonormal basis of L 2 (R), or, alternatively, a tight frame, as in (0) but with the ip^ k functions 
in place of ij)n,k- 

In principle, there are many ways (see below) of constructing wavelets (|^), but we will show in this paper that 
the method of quadrature mirror wavelet filters (QMF) has some features that set it apart from the alternative 
constructions. Several of the constructions are based on frequency subbands, and the subbands correspond to a 
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family of closed subspaces in the Hilbert space L 2 (R) , but we will show that this subspace structure is "optimal" for 
the QMF wavelets, in the sense that the subspaces cannot in a nontrivial way be refined into additional subbands. 
We will formulate this result in a mathematically precise fashion, which is based on a representation of the operators 
which define the QMF's. In fact, we give a formula for all these QMF's in the case of compactly supported wavelets. 

The present paper is concerned with the wavelet filters which enter into the construction of the functions 
tp( 1 ) ) tp( 2 \ . . . , t/i^ -1 ). These filters (see (0)-@ and ( |56| ) below) are really just a finite set of numbers which re- 
late the Z-translates of these functions to the corresponding scalings by x i— > Nx. Hence the analysis of the wavelets 
may be discretized via the filters, but the question arises whether or not the data which go into the wavelet filters are 
minimal. It turns out that representation theory is ideally suited to make the minimality question mathematically 
precise. (This is a QMF-multiresolution construction, and it is the minimality and efficiency of this construction 
which concern us here. While it is true, see, e.g., Refs. [|,[5|6|,|7],[| that there are other and different possible wavelet 
constructions, it is not yet clear how our present techniques might adapt to the alternative constructions, although 
the approach in Ref. || is also based on operator-theoretic considerations.) 

2. THE SCALING FUNCTION 

Since the wavelets come from the multiresolution functions, it is of interest to give explicit constructions for these. We 
do this in the present paper, which introduces two new tools for explicit constructions of multiresolution wavelet filters, 
(i) a certain infinite-dimensional loop group, and (ii) a certain family of irreducible representations of orthogonality 
relations (the Cuntz relations). Our viewpoint makes it clear, in particular, that compactly supported wavelets may 
be specified effectively with a finite set of spin vectors in C N (N is the wavelet scaling number, e.g., N = 2 in the 
dyadic case). Hence these wavelets are given by a finite set (in arbitrary configuration) of k Q-bits where it turns 
out that k is half of the length of the support of the wavelet in question. 

To better explain the minimality issue for multiresolution quadrature mirror (QMF) wavelet filters, we recall the 
scaling function ip of a resolution in L 2 (R) . 

Let g € N, and let do, ax, ■ ■ ■ , a 2g-i be given complex numbers such that 

2g-l 

E a ^ 2 ' ( ? ) 



fe=0 



and 



E 

k 



2 ifZ = 0, 

ak+2fak = 1 o if o. (8) 



In the formulation of (|^), and elsewhere, we adopt the convention that terms in a summation are defined to be zero 
when the index is not in the specified range: Hence, in (^), it is understood that a k +2i = whenever k and I are 
such thatJj + 21 is not in {0, 1, ... , 2g — 1}. When {a^ ; k = 0, 1, . . . , 2g — 1} is given subject to (0)-(||), then it is 
knowntmEj that there is a ip S L 2 (R) \ {0}, unique up to a constant multiple, such that 

23-1 

if(x) = E a k if(2x- k) , xeR, (9) 
k=a 

and if is of compact support; in fact, then 

supp (if) C [0, 2g - 1] . (10) 

(If H denotes the Hilbert transform of L 2 (M) , and ip solves ([)]) , then Hip does as well; but Hip will not be of compact 
support if ( |To| ) holds.) In finding if in (0), there are methods based on iteration, on random matrix products, and 
on the Fourier transform, see Refs. [l], ||, ||, and |Il| , |l2| , p"3 14; and the various methods intertwine in the analysis of ip, 



i.e., in deciding when ip (x) is continuous, or not, or if it is differentiable. 

Let ip be as in (||), and let Vo be the closed subspace in TL (:= L 2 (R)) spanned by {p (x — k) ; k <G Z}, i.e., by 
the integral translates of the scaling function ip. Let U (:= Un) be 

Uf(x):=N-*f(£), feL 2 (R), (11) 



Table 1. Discrete vs. continuous wavelets, i.e., £ 2 vs. I? (R) 
{0} < — • • ■ < — V 2 Vi V finer scales 
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the unitary scaling operator in H = L 2 (R). Then 

UV C V , (12) 

it is a proper subspace, and 

/\f/"Vo = {0}; (13) 

n 

see Ref. |9| and Ch. 5 of Ref. |l|. For N — 2, the situation is as in Table [l]. Setting 

V n := U n V , (14) 

and 



W„:=V n _ 1 GV„, (15) 

we arrive at the resolution 

V O = ]T W„, (16) 

7!>1 

and the wavelet function i\) is picked in Wo- We will set up an isomorphism between the resolution subspace Vo and 
I 2 (Z) , and associate operators in £ 2 (Z) with the wavelet operations in Vo C L 2 (R) . This is of practical significance 
given that the operators in £ 2 (Z) are those which are defined directly from the wavelet filters, and it is the digital 
filter operations which lend themselves to algorithms. 

For the general case of scale N (> 2) the space Vo UnVq splits up as a sum of orthogonal spaces W|*\ 
i = 1,2,..., N- 1. 



3. REPRESENTATIONS OF On AND TABLE | 
(DISCRETE VS. CONTINUOUS WAVELETS) 

The practical significance of the operator system in Table [t] (scale N = 2) is that the operators which generate wavelets 
in L 2 (R) become modeled by an associated system of operators in the sequence space £ 2 (:= £ 2 (Z) = L 2 (T)). (We 
will do the discussion here in Section ^| just for N — 2, but this is merely for simplicity. It easily generalizes to 
arbitrary N.) Then the algorithms are implemented in £ 2 by basic discrete operations, and only in the end are the 
results then "translated" back to the space L 2 (R). The space i 2 (R) is not amenable (in its own right) to discrete 



computations. This is made precise by virtue of the frame operator W : £ 2 (= L 2 (T)) — * Vo (c L 2 (M)) which may 
be defined as 

W:£ 2 3 (Z k )^Y,^(x-k)eL 2 (R). (17) 

fcez 

If the scaling function <p has then been constructed to have orthogonal translates, then W will be an isometry of £ 2 
onto Vo (c L 2 (M)). Even if the functions {<p (x — k)} keZ formed from p> by Z-translates only constitute a frame in 
Vo, then we will have the following estimates: 

/oo 
\(WO{x)\ 2 dx<c 2 Y,M\ (18) 
k^il fcez 

where c\ and C2 are positive constants depending only on the scaling function ip, or equivalently, 

c\ /2 -m\ P < nw^i i2 ( R ) < 4 /2 • iieiu . (19) 

Lemma 3.1. If the coefficients {ak ; k = 0, 1, . . . ,2g — 1} from ([)]) satisfy the conditions in (JsJ) , then the corresponding 
operator So: I 2 — > I 2 , given by 

( S o0k = ~^J2 ak - 21 ^ = ~?n H k £ Z, (20) 

p=k mod 2 

is isometric and satisfies the following intertwining identity: 

WS Q = UW, (21) 

where U is the dyadic scaling operator in L 2 (R) which we introduced in (j^. (Here we restrict attention to N = 2, 
but just for notational simplicity!) Setting 

b k := (-if a 2g -i- k , (22) 



and defining a second isometric operator Si : i 2 — > I 2 by formula (20) with the only modification that (bk) is used in 
place of (ak), we get 

S*S k =6 j<k l e * (23) 

and 

J2SjS] = lp, (24) 
3 

which are the Cuntz identities from operator theory^ and the operators Sq and Si satisfy the identities indicated 
in Table [l| 

Remark 3.2. For understanding the second line in Table [l], note that Sq is a shift as an isometry, in the sense of 
Ref. |l^, and C := S\£ 2 is a wandering subspace for So, in the sense that the spaces C, S C, S 2 C, . . . are mutually 
orthogonal in I 2 . To see this, note that ( p4| ) implies that 

(C :=) Sit 2 = £ 2 Q S £ 2 = kcr (S^) . (25) 

As a result, we get the following: 

Corollary 3.3. The projections onto the orthogonal subspaces in the second line of Table 1 which correspond to 
the Wi, W2, ■ • ■ subspaces of the first line (see (|l5|)) are as follows: 

proj C = SiSl = I - SqSq, 



on—In on— 1 o*n — 1 on Q*n 

P ro J *~>0 *- — D D — °0 D ■ 



Proof. Immediate from Lemma p.l| and Remark |3.2| . □ 

Remark 3.4. Any system of operators {Sj} satisfying (p3|)-(|2^) is said to be a representation of the C* -algebra O2, 
and there is a similar notion for On when N > 2, with On having generators Sq, Si, ... , Sn—i, but otherwise also 
satisfying the operator identities (p3|)-(p4]). 

Definition 3.5. A representation of On on the Hilbert space I 2 is said to be irreducible if there are no closed 
subspaces {0} Q Tio ^ I 2 which reduce the representation, i.e., which yield a representation of (^3|)-(^4|) on each of 
the two subspaces in the decomposition 



where 



£ 2 = H ® (I 2 e Ho) , 
£ 2 eH = (Ho) 1 - = {£ e I 2 ; (£ | 77) = 0, Vr? e H } 



(26) 



(27) 



Proof of Lemma 3J. Most of the details of the proof are contained in Refs. 17 and 0, so we only sketch points 
not already covered there. The essential step (for the present applications) is the formula (|2l| ) , which shows that W 
intertwines the isometry Sq with the restriction of the unitary operator U : f 1— > A=f (x/2) to the resolution subspace 

V C L 2 (M). We have: 



V pGZ VfcGZ / 



(by dl 
(by |)) 



V2- 



ip(x-p) 



(by 
(by I 



= (WSoO (x) 

for all £ € and all cc e E. This proves (||). 

The rest of the proof will be given in a form slightly more general than needed. For later use, we record 
the following table of operators on the respective Hilbert spaces L 2 (T) = I 2 and L 2 (R), and the corresponding 
transformation rules with respect to the operator W. Let N be the scale number, and let (cifc)j 
satisfying 



)k=0 



be given 
(28) 



and set mo (z) := Ylk=o a k zk 1 z € T. Then we have the following transformation rules. 
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TRANSLATION 

F(x) 1 — >F(x-l) 



£ 2 . £ 1 — > y^^k-Ni^i 

I Fourier transform ' 

L 2 (T): /^moW/(/) 



(60 
/(«) 



(6-1) 



real wavelets 
discrete model 



periodic model, 
T = M/2ttZ 



(29) 



Now the proof may be completed by use of the following sublemma, which we state just for N = 2. 
Sublemma 3.6. Let be a subspace of L 2 (T) which is invariant under multiplication by z 2 . Then there is an 
mi e L°° (T) such that & = {mi (z) f (z 2 ) ; / e L 2 (T)} and |E„;2=J m i H| 2 = l > a.a. z e T. 
Proof. The proof follows from the Bcurling-Lax-Halmos thcoremJl3 □ 



Completion of proof of Lemma 3T. Let mo(z) := -^^2k a k zk ■ Then we saw that Sof (z) := mo (z) / (z 2 ) is 
isometric in L 2 (T) , and the complementary space 



© := L 2 (T) G S L 2 (T) = ker (S* ) = /eL 2 (T); £ 

^ — Z 



mo M/ (w) = j 

then satisfies the condition in Sublemma |3.6j . Let mi be the function determined from the sublemma. Then, after 
multiplication by a suitable z 2p , we will get mi (z) — ^ fe b-^z k , with the coefficients 6^ as in (22). Setting 

S%f (z) := mi (z) f (z 2 ), it follows then that @-(||) are satisfied. □ 

Remark 3.7. The significance of irreducibility (when satisfied) is that the wavelet subbands which are indicated 
ill Table | are then the only subbands of the corresponding multiresolution. We will show that in fact irreducibility 
holds generically, but it does not hold, for example, for the Haar wavelets. In the simplest case, the Haar wavelet has 
g = 2 = N and the numbers 



a ° ai \ = I 1 1 
b h ll -1 



(30) 



Hence, for this representation of O2 on t 2 , we may take Ho = £ 2 (0, 1,2,...), and therefore Hq = f 2 (• ■ ■ , —3, —2, —1). 
Returning to the multiresolution diagram in Table [l], this means that we get additional subspaces of L? (M) , on top 
of the standard ones which are listed in Table |l|. Specifically, in addition to 

V n = U n V = WS£f (31) 

and 

W n = V n - 1 eV n = WS^- 1 S 1 £ 2 , (32) 

we get a new system with "twice as many", as follows: Vjf^ and W« , where 

V<+> = WS% (H ) , (33) 
W<+) = WS£- 1 S 1 (Ho) ; (34) 

and 

V<"> = WSS {K) , (35) 
Wt ] = WS*- 1 S 1 (H^) . (36) 

For the case of the Haar wavelet, see (pp|), 

V< +) cL 2 (0,oo), V^C^KO), 
or rather, Vo consists of finite linear combinations of ' ^-translates of 

\0 i/i€l\[0,l), y ' 

alias the step functions of step-size one, i.e., functions in L 2 (R) which are constant between n and n + 1 for all 
n G Z; and 

Vo +) =V ni 2 (0, 00) , V (_) =V nL 2 (-00, 0) . (38) 



Hence we get two separate wavelets, but with translations built on {0,1,2,...} and {...,—3,-2,-1}. In view of 
the graphics of the cascade approximation to the scaling function (see Refs. [ll] and [l2]), it is perhaps surprising that 
other wavelets (different from the Haar wavelets) do not have the corresponding additional "positive vs. negative" 
splitting into subbands within the Hilbert space L 2 (R) . 

Remark 3.8. There are other dyadic Haar wavelets, in addition to (|37|). For example, let 

if 0<x <2k + l, 



Vk f x ) = J V2k+i -■>--- ~ ■- ■ -> (39) 
\0 ifx£R\[0,2k+l). 

Then it follows that there is a splitting ofVa into orthogonal subspaces which is analogous to (p8[), but it has many 
more subbands than the two, "positive vs. negative", which are given in (p8|), and which are special to the standard 
Haar wavelet (|37|). For details on these other Haar wavelets, and their decompositions, we refer the reader to 
Proposition 8.2 of Ref. Eq. The corresponding m-functions of (p9[) are 



m (z) = -j=(l + z 2k + 1 ), mi (z) = ±=(l-z 2k+1 ), zeT. (40) 

Hence, after adjusting the O2~representation T with a (rotation) V € U2 (C), we have 

T f (z) = f (z 2 ) , T 1 f(z) = z 2k+1 f(z 2 ), feL 2 (T)=£ 2 , (41) 



and, of course, the two new operators Tq,T\ will satisfy the 02-identities (|23|)-(24), and the corresponding represen- 
tation will have the same reducing subspaces as the one defined directly from mo and mi. The explicit decomposition 
of the multiresolution subspaces corresponding to ( |38| ) may be derived, via W in Table |], from the corresponding 
decomposition into sums of irreducibles for the Oi~reprcscntation on i 2 which corresponds to (|3^). This means that 
the corresponding ( p6| ) which is associated with (^9|) and ([i"l| ) has more than two terms in its subspace configuration. 

4. WAVELET FILTERS AND SUBBANDS 

The power and the usefulness of the multiresolution subband filters for the analysis of wavelets and their algorithms 
was first demonstrated forcefully in Refs. [H] and ^o]; see especially p. 140 of Ref. [l9| and p. 157 of Ref. poL where the 
Oat -relations (p3|)~(p4|) are identified, and analyzed in the case N — 2. Around the same time, A. CoherO identified 
and utilized the interplay between £ 2 and L 2 (R) which, as noted in Section^ above, is implied by the 0Ar-relations 
and their representations. But neither of those prior references takes up the construction of Ojv-representations in a 
systematic fashion. 

Of course the quadrature mirror filters (QMF's) have a long history in electrical engineering (speech coding 
problems), going back to long before they were used in wavelets, but the form in which we shall use them here is 
well articulated, for example, in Ref. Ell Some more of the history of and literature on wavelet filters is covered well 



in Refs. 22 and 



The operators corresponding to wavelet filters may be realized on either one of the two Hilbert spaces I 2 (Z) or 
L 2 (T), T = R/27rZ, and L 2 (T) defined from the normalized Haar measure fi on T. But, of course, i? 2 (Z) = L 2 (T) via 
the Fourier series 



/(z) = ]T^ fc , (42) 



h-ez 

and its inverse 



& = /(*)= / z- k f(z)d^(z). (43) 
For a given sequence clq, ax, ... , ajv s -i, consider the operator So in £ 2 (Z) given by 



£ — > Sot, 



(44) 



and 



Setting 



Ng-l 



m 



a k z 



k=0 



and 



(S f)(z)=mo(z)f(z N ), /£L 2 (T), 
we note that So and So are really two versions of the same operator, i.e., 

(Soff = S (f) 



e T, 



(45) 



(46) 



(47) 



(48) 



when / = (£j.) from (^). (The first one is the discrete model, and the second, the periodic model, referring to the 
diagram (29).) Hence, we shall simply use the same notation So in referring to this operator in either one of its 
incarnations. It is the (|45|) version which is used in algorithms, of course. 

Let if £ 1? (R) be the compactly supported scaling function solving 



Then define the operator 



by 



Ng-l 

f{x) = dkf(Nx — k). 

fc=0 



W: f(Z) 



wt = j2tk<p(x-k). 



(49) 



(50) 



(51) 



The conditions on the wavelet filter {a,k} from Section [l| and in (p8|), may now be restated in terms of mo(z) in (Q) 
as follows: 



JV-l 



m (ze l «") 



fc=0 



= N, 



and 



mo(l) = VN. 



(52) 



(53) 



It then follows from Lemma 3.1 that W in (51) maps ^ 2 (Z) onto the resolution subspace Vo(C L 2 (R)), and that 

U N W = WS (54) 

where 

U N f(x) = N-if(jj), / e L 2 (R), iel. (55) 
We showed in Ref. |§| that there are functions mi, . . . , m^_i such that the N-by-N complex matrix 

1 / -M,r \N-1 

(m 3 (e l — z)\ (56) 



is unitary for all zeT. If wc define 

S j f(z)=m j (z)f(z N ), feL 2 (T), zeT, (57) 



then 



and 



S*Sk — <5 3 ,fe/ L 2( T ), (58) 



N-l 



£ SjS* = I LHT) . (59) 



(Condition ( |53|) is not needed for this, but only for the algorithmic operations of the Appendices in Refs. [ll] and 
|l2| .) We also showed that the solutions (mj)j^ 1 to (|56| ) are in 1-1 correspondence with the group of all polynomial 
functions 

A: T — >Uat(C) (60) 
where U/v(C) denotes the unitary N x N matrices. The correspondence is m <-> A with 

N-l 

m j {z) = Y J A j , h {z N )z k , (61) 

fc=0 

and 

A jt k(z) = ^ E w- k mj {w), (62) 

W N =Z 

This correspondence plays a central role in the proofs in Refs. 11 and [l2| Wc also show in Ref. || that if m,Q is given, 
and if it satisfies (|52]), then it is possible to construct m 1; . . . , tojv-i such that the extended system m , m 1; . . . , m^r-i 
will satisfy ([56]). As a consequence, A in (|6|) will be a U/v(C)-loop, i.e., A: T — > Ujv(C), and moreover, the original 
mo is then recovered from (^lj) for j ' = 0. 

By virtue of (^8|)- (^9|), L 2 (T), or equivalently ^ 2 (Z), splits up as an orthogonal sum 

SjifiZ)), j = 0,1,..., N -1. (63) 

We saw that the wavelet transform W of (pp|)-(pT|) maps £ 2 (Z) onto Vo, and from (54) we conclude that W maps 
Sq(£ 2 (Z)) onto Un(Vq) (=: Vi). Hence, in the TV-scale wavelet case, transforms the spaces Sj(l 2 (Z)) (c ^ 2 (Z)) 
onto orthogonal subspaces , j = 1, . . . , JV — 1 in L 2 (R), and 

JV-l 

Wi = v e Vi = V ® Wp , (64) 



where 



Wj j) =^£ 2 , j = l,...,iV-l. (65) 



Each of the spaces Vi and W{ is split further into orthogonal subspaces corresponding to iteration of the operators 
So, Si, ... , Sn-i of (p^)-(p^)- It is the system {Sj}^^ 1 which is called a wavelet representation, and it follows that the 
wavelet decomposition may be recovered from the representation. Moreover, the variety of all wavelet representations 
is in 1-1 correspondence with the group of polynomial functions A as given in (|60|). The correspondence is fixed 
by (|6l|)-(62). Operators {Sj} satisfying (]58|)-(p9|) are said to constitute a representation of the C*-algebra On, the 
Cuntz algebrajl3 and it is the irreducibility of the representations from (|57j) which will concern us. If a representation 



([57]) is reducible (Definition 3.5), then there is a subspace 

0C Wo C£2 (T) (66) 



which is invariant under all the operators Sj and S*, and so the data going into the wavelet filter system {nij} are 
then not minimal. 



5. THE MAIN THEOREM 



The main result will be stated in the present section, but without proof. Instead the reader is referred to Ref. [T^ for 
the full proof, and for a detailed discussion of its implications. We noted above that the representation ( p7\ ) given 

(A) 

from a QMF system rrij = rrij 

not admit further reduction into a refined system of closed subspaces of L 2 (K) 
Theorem 5.1. The representation 



in , ' via (^TJ) — (|6^) is irreducible if and only if the subbands are optimal, in that they do 



8) A) f (z) = mf > (z) f (z N ) , / e L 2 (T) , z e T 

is an irreducible representation of On on L 2 (T) if and only if A: T — > ~Un (C) does not admit a matrix corner of 
the form 



V 



(i 











(67) 



for some V € Um (C), and where no, Tlx, ■ ■ ■ , tim-1 G {0, 1,2,...}. 
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